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ABSTRACT 


The solution of optimum rendezvous trajecto- 
ries with bounded thrust magnitude and interim 
coasting arcs has been obtained. The Lagrange 
multipliers, determined from the impulsive solution, 
are used as starting values. The three-dimensional 
equations have been derived for the assumptions of 
a mass particle moving in an inverse square force 
field in a vacuum. A practical solution time was ob- 
tained by the use of a convergence technique which 
employs the impulsive solution to generate correc- 
tions to the initial multipliers. 


ii 



THE UTILIZATION OF AN IMPULSIVE SOLUTION IN 


SOLVING OPTIMUM TRAJECTORY PROBLEMS 

WITH BOUNDED THRUST MAGNITUDE 

By Donald J. Jezewski 
Manned Spacecraft Center 


SUMMARY 


A method is presented for obtaining solutions to nonlinear optimal trajectory 
problems with bounded thrust magnitude by use of the corresponding impulsive solution 
to generate the starting values of the Lagrange multipliers. The equations were de- 
rived in three dimensions for a mass particle moving in an inverse square force field 
in a vacuum. The convergence scheme which makes use of the impulsive multipliers 
has an efficient solution time as compared to a first-order perturbation technique. 

An example problem of a transfer from Earth to Mars for fixed initial thrust ac- 
celerations has been used to demonstrate the solution technique. The number of iter- 
ations and the time required for a solution increased rapidly as the assumptions of the 
impulsive solution were violated. 


INTRODUCTION 


The solution of optimal trajectories with bounded control variables, which uses 
the indirect method of the calculus of variations, has been treated by many investiga- 
tors (refs. 1 to 5). This formulation results in a two- point boundary -value problem 
which usually has unknown initial or final values of the Lagrange multipliers. The pro- 
cedure has simply been to make an intuitive guess of these values. For special prob- 
lems in particular coordinate frames, an experienced investigator can attach some 
physical meaning to these multipliers (ref. 6). The extreme sensitivity of the terminal 
boundary conditions to variations in the initial or final values of these variables is well 
known. This deterrent to the indirect solution method gives rise to other techniques 
such as the gradient or steepest- ascent method (ref. 7). Difficulties are also inherent 
in this method; for example, the choice of proper step size and the manner in which it 
is varied (ref. 8). The tendency of investigators is to look to other methods of solu- 
tion. One interesting possibility, a hybrid optimization technique, consists of the 
method of steepest ascent for the initial phase and one of the indirect methods for the 
latter phase. 



Pines (ref. 9) suggested that the Lagrange multipliers obtained by assuming an 
impulsive solution be used as starting values. He reasoned that, "We can look at the 
impulsive solution as a limiting point in a simply connected region in the space of the 
initial conditions of the adjoint variables. " This indicates that, at worst, a solution 
could be obtained by stepping the initial thrust acceleration from the impulsive case to 
the desired finite level. 

Handelsman (ref. 4) solved two-dimensional Earth-to-Mars trajectory problems 
by use of the impulsive multipliers as starting values. The author's work is an ex- 
tension of and improvement on Handelsman' s solution. The extension involves adding 
the third dimension to the solution. The improvement derived is the increased effi- 
ciency, as defined by the number and time for an iteration, in obtaining a solution. 
This efficiency is accomplished not only by using the impulsive multipliers as starting 
values, but also in the convergence scheme of a finite thrust-acceleration problem. 
The method of using an approximate closed-form solution in a convergence scheme 
(feedback loop) of a nonlinear problem is exemplified in Baffin's work (ref. 10) and in 
Jezewski's work (ref. 6) for the circumlunar and time-optimal trajectory problems. 


SYMBOLS 


A.. 

i] 


A 


z 


B 

C l> C 2’ ' • • ’ C 6 


e 

f f f 

l r 2 ’ • ■ • ’ 6 



k 


m 

P 


q 

r 


matrix defined in equation (78) 

azimuth angle measured from 1^ , deg 

terminal boundary vector 
constants defined in equation (78) 

constant defined in equation (19) 

eccentricity of the coasting ellipse 
functions defined in equations (51) to (56) 

specific impulse, sec 

logical control switch 
vehicle mass, slugs 

period of circular orbit whose radius is r^, sec 
function defined in equation (57) 

radius to mass particle from center of attracting body, ft 


2 


transformation matrix defined in equation (59) 
thrust, lb 
time, sec 

components of total velocity in directions in 1^, 1^, and 

rections, respectively, ft/sec 

total velocity, ft/sec 

exhaust velocity, ft/sec 

nondimensional exhaust velocity 

inertial reference frame, ft 
nondimensional velocity 
unit vector in the v direction 

unit vector in the direction normal to the plane of motion 

unit vector in p, <fi, and 6 directions, respectively 

included angle between thrust and velocity vectors, deg 
variable defined by equation (103) 
flightpath angle, deg 
change of 

nondimensional time rate of change of latitude 
included angle between 1^ and 1^ unit vectors, deg 

latitude angle, deg 

switch function defined in equation (17) 
primer vector 
Lagrange multipliers 


mass fraction 



v 

TT 

P 

O 

T 

0 

X 

V 

n 


3 / 2 

gravitational constant, ft /sec 

true anomaly angle, deg 
3. 1415927, rad 
nondimensional radius 

nondimensional time rate of change of longitude 
nondimensional time 
longitude angle, deg 

thrust pitch angle measured from local horizontal plane, deg 
thrust yaw angle measured from 1^, deg 

time rate of change of e, sec * 


Subscripts: 

C 

I 

1, j 

max 

R 

U, V, w 

0,1 

2, f 
P 


refers to a corrected solution 
refers to an integrated solution 
indices 

maximum value of 
reference quantity 

refer to components in 1^, 1 , and 1^ directions, respectively 

refer to conditions before and after first impulse 
refer to conditions before and after second impulse 

in the direction of the unit vector 1 

P 

below a variable indicates a vector quantity 
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Superscripts: 




a, b 

value of a quantity after 

and before an event 

R 

with respect to the 1^, 

1 R , and 

1 ,p coordinate frame 

S 

with respect to the 1^, 

V and 

1q coordinate frame 

-1 

inverse of a matrix 



A 

approximate value of variables 



Operators: 

( ) derivative with respect to time t 

( )' derivative with respect to time r 

OPTIMUM TRAJECTORY EQUATIONS 

The nonlinear, nondimensional equations of motion of a mass particle moving in 
an inverse square force field in a vacuum and acted upon by a thrust acceleration can 
be written in the following manner (see appendix): 

Z y ’ - p(a 2 cos 2 0 + 6 2 ) + + jj v 0 sin x = 0 (1) 

n ' 


per’ cos 9 

+ 2 a (z cos 9 - p6 sin 9 j + — v 
l v H ) p e 

cos x cos = 0 

(2) 

p5' + 

2 n< 

26 Z + pa sin 9 cos 9 + - v cos 
v n e 

X sin i// = 0 

(3) 


1 

N 

< 

II 

O 


(4) 


Q 

II 

O 


(5) 


0’ - 6 = 0 


(6) 


5 


p' +(3 = 0 


(7) 


The prime refers to the nondimensional time t = tV-^/r and (3 is proportional to the 

thrust T. The state variables are the position and velocity coordinates p, <P, 6, 

Z , a, 5, and the mass p. The control variables are the magnitude and direction of 

the thrust (3v , x> and i//, respectively. The exhaust velocity v is assumed con- 

stant. The axis system and the associated notation is illustrated in figure 1. 




Figure 1. - Coordinate system and angle definition. 


For this system of constraint equations, it is desirable to determine the thrust-vector 
program (bounded in magnitude) for a trajectory between two fixed positions, veloci- 
ties, and times which minimizes the propellant consumption. The necessary and suf- 
ficient conditions for an optimum trajectory may be found in reference 1. It will suffice 
to state the variational equations and indicate their functions. 


A^' = 2 (^ 2 ° cos ^ + ^3^) " ^4 


(8) 


V = _2ct ( x i cos 0 - x 3 sin e ) + x 2 hr - 6 tan 7 ' 


(9) 


Z A 

A 3* = " 2 (V + V Sin 0 ) + X 3 ' ~p 


( 10 ) 


6 



. /c2 2 2 0 2 

= -X^/S + cr cos 0 + 


+ Xg^'cos 0 - 2cr6 sin 0) 


+ X^S' + a^sin 0 cos 0) 


v = ° 


( 11 ) 

(12) 


2 2 
\q = 2Xjpa sin 0 cos 0 + X^pa cos 20 

-Xgjpa'sin 0 + 2a sin 0 + p5 cos 0^J 


(13) 


V - HW 


(14) 


Xj Xg 

tan x = t~ cos tan \p = t— 


(15) 




(16) 


K = 7T W - x , 




(17) 


0=< 


r /3 , k> 0 

max’ 

0 , k< 0 


(18) 


C = -X n Z ’ - X 9 pa’ cos 0 - X„p6* - X.p' 


- V - V - V' 


(19) 


Equations (8) to (14), the Euler- Lagrange equations, are the first necessary conditions 
for an optimum trajectory. The optimum conditions on the thrust vector are given by 
equations (15) to (18). These equations show the thrust direction to be alined with the 
primer vector A (Lawden's notation, ref. 2) and indicate the condition for the change 
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in magnitude of the bounded thrust (the switching function k (t)). The condition for 
k 5 0 (intermediate thrust arcs) will not be considered in this study. The primer vec- 
tor A has components X^, Xg, and Xg in the directions of the unit vectors 1^, 1^, 

and 1^ , respectively. The direction cosines of the thrust vector can be written as 


smx= m 


cos X COS 1// = pq- 


cos X sin # = pq- 


( 20 ) 


( 21 ) 


( 22 ) 


This allows a replacement of the angles x and in the formulation by the primer- 
vector magnitude and its components. Equation (19) results from a first integral; C Q 

is a constant if the external accelerations are explicitly independent of the variable r . 


IMPULSIVE RELATIONSHIPS 


An impulsive solution is assumed to occur in two possible manners: the external 
force acceleration approaches infinity or the independent variable t goes to zero. 

The form of equations (1) to (14) is unsuitable for determination of the relationships be- 
tween the variables when an impulse is applied. A change in variables is made from 
r to p because the mass ratio has a finite variation across the impulse. Performing 
this change in the independent variables and taking the limit of equations (1) to (14) as 
/3 approaches infinity, the following equations result. 


AZ y = |AV| sin X 

(23) 

p Act cos 9 = |AV| cos x cos j// 

(24) 

p A6 = |AV| cos x sin 

(25) 

Ap = 0 

(26) 

A<p = 0 

(27) 


8 



& 

i! 

O 

(28) 

AA, = 0 

(29) 

aa 2 = o 

(30) 

AAj.O 

(31) 

A A. = A. Act cos 9 + Ag A6 

(32) 

aa 5 = o 

(33) 

AXg = -pX 2 Acr si n $ 

(34) 

A V V eW(;^) 

(35) 


The quantity AV is 



Jh 

P 


This equation is not actually used to determine the magnitude of the impulse because 
the change in the mass ratios is not known a priori. The AV is determined by dif- 
ferencing the velocity vectors before and after the impulse. 


a Y=V 1 -V q (37) 


Because the impulse is alined with this vector, the thrust direction angles x and \}/ 
are also specified. 
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In equations (23) to (35), the quantities p, <p, 6, A^, Ag, and Ag do not 

change, which indicates that the position and primer vector are invariant over the im- 
pulse. The variable Ag is a constant of motion of the solution and is unaffected by the 

impulse. The relationships between the variables that change across the impulse are 


z a = z b + |av| Sin X 


a b M 

a = a + „ cos x cos xu 

p cos 6 A ^ 


6=6 + cos v sin xLf 

p A 


b l A -l l A -l 

A^ a = A^° + Ag cos x cos rp + A g cos x sin ip 


A^ a = Ag b - AgjAVj cos X cos ^ tan 0 


A a = A b 


\M I* ) 


where equations (39) and (40) were used in equations (41) and (42). 


PRIMER VECTOR ON A COASTING ARC 


Lawden's equations for the Lagrange multipliers and their time rate of change 
on a coasting arc (ref. 2) are 


Aj = Cj cos v + Cg e sin v + C^M 


A 2 = CjfgM + CgfjM + C 3 f 3 (^) + 
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(46) 


where 


and 


X 3 “ 


^Cg cos v + Cg sin vj 




Xj = v j^-C^ sin ^ + Cg e cos v + 


) n 

“e + cos v 

r 

s + COS V 

fl (V) 


- Cg e sin v 


+ Cr 


COS V . , s 

rrr. - h M 

fj O') 


e sin v 

i*(y) 


v [-C 5 sin v + Cg(e + cos y)J 


X 3 “ 




. ( L 

v - (^ct cos 9 


2 2 r 


+ s 2 ) 


1/2 


fj(i') = 1 + e cos v 


ef„(^)sin v 
-cos v 6 V ' 


r / \ V 


, , , f i M v v) 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


(53) 
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t / \ ■ ef e (y)cos v 

f A (v) = sin v 6 V ' 

■ “RT 


(54) 


- - [* + qW] sin " 




sin v (l + 2e^)sin v 
(l-e 2 )fjM 

fie tan" * [q tan ($)] 

' G-e 2 )^ 


( 55 ) 


( 56 ) 


_ f(l - e)l 

- LT^r^rJ 


1/2 


( 57 ) 


This set of equations is a function of only three quantities: the constants 

Cj, Cg, . . . , Cg defined by the initial conditions on the multipliers and their deriva- 
tives, the true anomaly v, and the eccentricity e of the coasting arc. The compo- 
nents Ag, and of the primer vector as defined above are with respect to an 

altitude, range, and track coordinate frame. The transformation to the spherical co- 
ordinate frame is given as 


A 


S 



( 58 ) 


where the superscripts R and S refer to the altitude, range, and track coordinate 
frame, and the spherical coordinate frame, respectively. The transformation S is 
given as 


S = 


10 0 
0 cos e -sin e 

0 sin e cos e 


(59) 


where e (fig. 2) is a negative rotation about the unit p- axis. 
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z 



Figure 2. - Angular relationship between an altitude, range, and track frame 

and a spherical coordinate frame. 


The time rate of change of this vector is given as 


; S • R nu .S 
A — SA + X A 


(60) 


where the vector O is the angular rate between the two coordinate frames. This vec- 
tor is directed along the p-axis and may be readily verified to have a magnitude given 
as 


|fi| .= = -cr sin 0 (61) 

INITIAL MULTIPLIER DETERMINATION 


It may be noted from equations (29) to (31) that the components 
of the primer vector are continuous across an impulse. Therefore, 


X 2 , and X 3 ) 


X x b = X x a (62) 

X 2 b = X 2 a (63) 
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( 64 ) 




a 


From the previous section the primer vector and its derivative on the coasting arc are 
known by equations (44) to (49). This vector is, therefore, also known before the im- 
pulse. 


The Euler- Lagrange equations are a homogeneous set. This condition indicates 
that only their ratios are significant and that one of the multipliers may be chosen ar- 
bitrarily. Since the multiplier X^ will remain nearly constant over the solution, its 

initial value will be chosen on the following basis: 



-}<«<} 




(65) 


where 


-1 

a = cos 


T • V 

nrra 


( 66 ) 


This choice of sign for the initial value of X^ is based on whether the impulse is a 
posigrade or retrograde maneuver. Using equations (9) and (33), the constant Xg 
may be determined as 

. b s a a 

Xg = Xg = p cos 6 

-2ct^Xj cos 6 


X„l 6 tan 6 

2 \P ) 


- Ag sin 


V 


(67) 


Similarly, from equations (10) and (34) 


Xg S = XgZ v - 2p^XjS + XgCT sin 9^j - pXg' 


( 68 ) 


The initial value of the mass 


multipliers X 




must be determined before a solution 
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can be generated. Recalling equation (19) 


C = -X..Z ' - X-po’cos 9 - X„p6’ - X p’ 

O IV ^ o *± 

- v - V - V 

and using the equations of motion, this constant may be written as 


C o = X l 


-p(cr"cos“0 


2 2 . 


.2^ 1 

+ 6 / + ~2 


+ 2X2^^Z v cos 9 - pd sin d'j 


X 0 (26Z + pcr^sin 9 cos 0) -X.Z - X c cr - X„6 
3\v K / 4 v 5 6 


fix 


(69) 


where equation (17), the relationship for k, has also been used. During the coasting 
portion of the solution, k < 0 and fi = 0, which facilitates a solution for the constant 
C in terms of known functions evaluated after the impulse. Since C has now been 

uniquely defined, the initial value of k , and hence of X , may be determined from an 

M 1 

evaluation of equation (69) before the impulse. Solving for X 


x b * 

p /3 


f 1 

C o + n lAl +X i 

(2 2 n .2^ 1 

p\a cos 0 4-oj 2 


P _ 


- 20 X 2 ^Z v cos 9 - p6 sin 0^ - X g ^26 Z^ + pcr^sin 0 cos 0^ 


4- X . Z 4- X_CT 4* X 6 
4 v 5 6 


(70) 


The initial values of the Lagrange multipliers have now been determined as func- 
tions of the state variables before and after the impulse, of the magnitude of the im- 
pulse, and of the primer vector and its derivative after the impulse. The determination 
of this vector and its derivative requires a solution of the constants C^, Cg, . . . , C g 

in Lawden's equations. A solution for these constants requires six independent 


15 



equations which govern the relationship of the multipliers on the coasting arc. The 
simplest set of equations available is 


cos - X 20 tan x 0 = o 

(71) 

cos i// { - X 2f tan x f = 0 

(72) 

X 30 “ X 20 tan ^0 = 0 

(73) 

X 3f ” X 2f tan = 0 

(74) 


X 10 X lf 


sm X Q sin x f 


(75) 


f |av| ^ 

a ± 1 = X 20 I 2ctcos 9 + — — cos x cos 


5 b M 


+ X 30\ 26 + COS X Sin 


(76) 


where the following notation has been used 


^io-^i^o) 


A ^ — A 


x if _ x i( t f) 


(77a) 

(77b) 


Equations (71) to (74) result from evaluation of equation (15) at the initial and terminal 
values of the coasting arc. Equation (75) was obtained by evaluating equation (20) at 
both impulsive times and by eliminating the magnitude of the primer vector between 
the two relationships. This is possible because the magnitude of the primer vector 
must take on identical values at the two impulsive times, as is indicated by examina- 
tion of equations (14) and (17). Equation (76) is equation (8) evaluated after the first 
impulse, with having been eliminated by use of equations (41), (65), and (66). 
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After transformation of this set of equations into the altitude, range, and track 
coordinate frame by means of equations (58) to (61), the constants C^, Cg, . . . , Cg 

can then be evaluated by inversion of the matrix of coefficients A. ^ and by solution of 

the following equation: 


c i - A i, 


-1 


0 

0 

0 

0 

0 

±1 


i, j — 1, 2, . . . , 6 


(78) 


The proper sign is determined from the relationships established in equations (65) and 
(66). It is noted that the values of the A^ and Ag multipliers after the first impulse 

were required for a solution of the initial values of the Lagrange multipliers. 


CONVERGENCE TECHNIQUE 


With a technique for generating the initial values of the Lagrange multipliers 
from an approximate closed-form solution, a convergence process may be established 
for the exact solution. The basis for convergence is that the impulsive (approximate) 
and nonimpulsive (exact) trajectories have solution curves which are similar, such that 
an identical change in the initial values of each will produce approximately the same 
change in the terminal boundary conditions. These terminal boundary errors can be 
nullified by generation of corrections to the initial values from the impulsive solution. 
Convergence will occur if assumptions of the impulsive solution are not violated ex- 
cessively. 

Figure 3 is a logical flow chart of the iterative technique used in the convergence 
process. A switch k governs the various phases of the solution. 

With the desired terminal boundary conditions and the switch k equal to 0, the 
impulsive solution is solved for the initial values of the Lagrange multipliers Ar and 

is denoted as the reference solution. If the nonimpulsive system of equations is inte- 
grated to a terminal cutoff using the initial values A , the resulting terminal vec- 

— XV 

tor Bj should resemble the desired terminal boundary conditions. The error in these 

conditions is indicative of the accuracy of the approximate solution. After replacement 
of the desired boundary conditions by those obtained from the nonimpulsive solution 
B = Bj the impulsive solution is resolved for the initial vector A. The switch k is 

set equal to +1. In the impulsive solution, a change in the initial vector of 


AA = A - A r (79) 
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Bj Resulting integrated vector 

X Lagrange multiplier 

Ar Reference Lagrange multiplier 

Ac Corrected Lagrange multiplier 


Figure 3. - The use of an impulsive solution in a convergence technique, 
has produced a change in the terminal conditions of 


AB = Bj - B 


(80) 


Because the impulsive solution is an approximation to the exact system of equations, it 
is assumed that the negative of this correction should produce an approximate change 
in the terminal boundary conditions of the nonimpulsive solution. Therefore, the cor- 
rected initial values of the Lagrange multipliers are 


— C 




(81) 


The switch k for the next and all future passes is set equal to -1. The procedure is 

repeated with one exception. The correction to the reference vector is made with 

— K 
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respect to the previous integrated solution X„. This correction is given by 


AX — X — Xq 


( 82 ) 


The solution is continued in this manner until the terminal boundary conditions ap- 
proach the desired conditions within some arbitrary tolerance. 


EXAMPLE PROBLEM 


To adequately demonstrate this solution technique, a fuel-optimum Earth-to-Mars 
transfer problem is solved for a wide range of initial thrust accelerations. The solu- 
tion is initiated on Julian date 2 442 680 with a 213-day transfer time. The exhaust ve- 
locity is fixed at 1. 0231545 EMOS (Earth mean orbital speed). Convergence was rapidly 
obtained for thrust accelerations ranging from the impulsive case to a value of 

5. 0 x 10 ^g (lg = 32. 2 ft/sec^). Solutions for accelerations less than this value 
could not be obtained without modifying the covergence technique. The variation in the 
initial values of the Lagrange multipliers with initial thrust acceleration for this trans- 
fer problem is illustrated in figure 4. The multipliers can be noted to diverge from 
their impulsive values as the initial thrust acceleration is decreased. It may appear 
that the initial values of the Lagrange multipliers are insensitive to the initial thrust 
acceleration. In this respect, this figure is misleading. The sensitivity of the mul- 
tipliers to changes in the initial conditions is obscured by the nondimensionalization 
of the problem variables. Because, for this solution technique, only the mass mul- 
tiplier X is a function of T/m and the remaining lambdas use the impulsive solution 

multipliers as starting values, convergence will be limited to thrust-acceleration levels 
which do not violate excessively the assumptions of the analysis. For example, the 
impulsive solution is based on a thrusting logic cycle of thrust-coast-thrust. If, for a 
finite acceleration level, the starting values produce a solution which does not complete 
this cycle before reaching a terminal cutoff condition, convergence will not occur. 

The significance of figure 4 should be reemphasized concerning the problems to 
which this solution technique is applicable. It will be recalled that the equations de- 
veloped in this analysis were formulated in a nondimensional form. In this formula- 
tion, the particular example problem is immaterial with respect to the developed 
equations. The characteristics of the Lagrange multipliers for the investigated range 
of the initial thrust acceleration are shown in figure 4. 

A measure of the efficiency of the solution technique is the number of iterations 
and the time required for a convergence to a desired set of boundary conditions. An 
iteration is defined as the number of corrections that are made to the multipliers such 
that the terminal constraints are satisfied on the succeeding solution within an arbi- 

-4 

trary tolerance. For this example problem, a uniform tolerance of 2. 0 x 10 was 
selected because all variables are nondimensionalized. The data for each iteration 

-4 

of the Earth-to-Mars transfer with an initial thrust acceleration of 9. 0 x 10 g are 
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Initial thrust acceleration/ g, nondimensional 



TABLE I. - EARTH- TO- MARS TRANSFER (213-DAY) a 





Iteration no. 




V <iria.UJ.eo 

0 

1 

2 

3 

4 

5 

6 

V p , EMOS .... 

0. 70447625 

0. 70764348 

0. 73613807 

0. 73759156 

0. 7394565 

0. 7394178 

0. 7395386 

r F , deg 

3. 5489735 

2. 6832596 

1.4632210 

1.3113635 

1. 236200 

1. 233388 

1. 2288496 

A zf , deg ... . 

90. 521965 

90. 098216 

90. 118342 

90. 123595 

90. 12412 

90. 12409 

90. 12396 

Tp,, AU • • • • 

1. 7477922 

1. 7384044 

1. 6700365 

1. 6666388 

1.662387 

1.6624722 

1.6621979 

0 F> deg 

2. 5647125 

1. 8270476 

1. 9067470 

1. 8395678 

1. 849077 

1.844119 

1.8451757 

deg 

143.07035 

142. 06986 

143. 11270 

143. 15160 

143.2138 

143. 2086 

143. 21273 

% ' 10 * 

2. 8987639 

3. 2015208 

3. 1592476 

3. 1478997 

3. 145361 

3. 144586 

3. 144495 

x 20 10 • ■ • • 

5. 2501955 

5. 1380339 

5. 1221021 

5. 1420125 

5. 1414238 

5. 1427927 

5. 142652 

x 30 • 10 1 .... 

2. 5655474 

1. 7708472 

2.0219506 

1. 9684969 

1.990323 

1.985170 

1. 987014 

»50 ' 101 

. 24761119 

-.09333626 

. 00145934 

.00513479 

.01298817 

. 01273934 

.01327756 

x 60 • 10 1 .... 

-5. 3477658 

-3.6911911 

-4. 2379409 

-4.092365 

-4. 1393818 

-4. 126159 

-4. 130388 

Xpi • 10 1 .... 

0 

6. 4003868 

6.0136106 

6. 1206711 

6. 1079195 

6. 117931 

6. 1162759 

6. 1170768 

t, days 

219.61825 

219. 15351 

213.76202 

213.37727 

213.02954 

213. 02164 

213. 00003 


a Julian date 2 442 680. 



listed in table I. It is noted that convergence to the desired boundary conditions is 
rapid and uniform. Computer time for this solution on a CDC 3600 electronic data 
processing machine was approximately 2 seconds. 

-4 

An initial thrust acceleration of 5. 0 x 10 g required 20 iterations and a solution 
time of approximately 7 seconds. This indicates that, as the assumptions of the impul- 
sive solution are violated, the number of iterations and the time required for a solution 
rapidly increase. 

For a comparison with other convergence techniques, an example problem with 

an initial thrust acceleration of 3. 1 x 10 g was solved using a first-order perturbation 
technique. This convergence method required six times as much computer time as the 
present technique. 


CONCLUSIONS 


A practical technique has been presented for obtaining solutions to optimum tra- 
jectory problems with bounded thrust magnitude by use of an impulsive solution to 
generate the starting values of the Lagrange multipliers. 

This simplifies the problem by eliminating the necessity of choosing the initial 
values of the Lagrange multipliers. The convergence scheme employs the impulsive 
solution to generate corrections to the initial multipliers. The technique is demon- 
strated by solving an Earth-to-Mars transfer problem for a wide range of initial thrust 
accelerations. Computer (CDC 3600) solution time for the specific case of the thrust- 

-4 

to-mass ratio = 9. 0 x 10 g was 2 seconds. The solution time was shown to 

increase rapidly as the assumptions of the impulsive solution are violated. In a com- 
parison case with a T/m o = 3. 1 x 10 ^g, the convergence scheme was six times more 
efficient than a first-order perturbation technique. 

Lawden's equations for the Lagrange multipliers were used on all nonthrusting 
arcs. This facilitated a rapid solution of the Euler- Lagrange equations and eliminated 
the need for integrating on coasting arcs. 

For a problem in which the variables have been nondimensionalized, the initial 
value of the Lagrange multipliers appear to be constant for a wide range of value of 

7 m o- 
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APPENDIX 

DERIVATION OF NONDIMENSIONAL EQUATIONS 

The equations of motion in spherical coordinates of a mass particle acted upon by 
the accelerations of gravity and thrust T/m are 


. (u 2 + w 2 ) M G T . 

v - + — sy- - — sin x = 0 

r 2 m 

r 

(83) 

u(v - w tan 0) T , n 

u + — 1 cos x cos xb = 0 

r m ^ 

(84) 

( 2 ^ 

\vw + u tan 6J T , n 

w + -X l cos y sin ib = 0 

r m a ^ 

(85) 

r - v = 0 

(86) 

r <p cos 6 - u = 0 

(87) 

Qi. 

1 

II 

O 

(88) 

T 

m +Y -=0 

(89) 


e 


where x and ip are arbitrarily chosen control variables and V g is the exhaust ve- 
locity, assumed constant. Choosing a reference position r and a reference velocity 

K 

V.^, defined as circular velocity at the position r^, 


1/2 


V. 


R 


_G 

: Ry 


(90) 
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the nondimensional equations of motion are 


(z 2 + z 2 ) 1 , 

P 

(91) 

fz - Z tan 0) f 

Z u +Z u p “ + p v e cos x cos 0=0 

(92) 

(z Z + Z 2 tan 0) , 

„ . \ v w u / U T , 

Z w + p + j;v e cosxsin^-0 

(93) 

p’ - z v = 0 

(94) 

P0 T cos 0 - Z =0 
u 

(95) 

p0‘ - Z =0 
w 

(96) 

+ 

II 

O 

(97) 

The superscript prime refers to the nondimensional time 


2* V 

T - p _ r 

R 

(98) 

where P is the period of the circular orbit. The new variables are defined by the 
following equations. 

il 

(99) 

r 

P = — 

R 

(100) 
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(101a) 


Z = ^ 
v V. 


R 


Z = ^ 
u V. 


R 


(101b) 


Z = JL 
w 


(101c) 


P = 


m 

m 


( 102 ) 


0 = - 


m r R 
m o V R 


(103) 


The equations of motion are further transformed by letting 


a = <p' = 


u 


p cos 0 


(104) 


5 = 0' = 


w 


(105) 


This change of variables results in the following nondimensional equations in terms of 
the variables Z y , a, 5, p, 4>, 0, and p. 


Z ' - p(cr 2 cos^0 + 6 2 ) + At + 77 V sin X = 0 

V Z U 6 

p 


(106) 


per' cos 0 + 2a ^Z v cos 0 - p6 sin 0^ + ^ v e cos X cos ^ = 0 


(107) 


2 u* 

p6' + 26Z + pa sin 0 cos 0 + — v cos x sin ^ = 0 
v fi e 


(108) 
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p’ - Z =0 

V 

(109) 

1 

Q 

II 

O 

(110) 

0’ - 6 = 0 

(HD 

H’ + j8 = 0 

(112) 
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